Multiband effects on the conductivity for a multiband Hubbard model 
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The newly discovered iron-based superconductors have attracted lots of interests, and the corresponding the- 
oretical studies suggest that the system should have six bands. In this paper, we study the multiband effects 
on the conductivity based on the exact solutions of one-dimensional two-band Hubbard model. We find that 
the orbital degree of freedom might enhance the critical value C,. of on-site interaction of the transition from 
a metal to an insulator This observation is helpful to understand why undoped High-rj. superconductors are 
usually insulators, while recently discovered iron-based superconductors are metal. Our results imply that the 
orbital degree of freedom in the latter cases might play an essential role. 

PACS numbers: 71.30. +h, 71.10.Fd, 74.70.-b 
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Recently discovered iron-based superconductors have at- 
tracted lots of experimental and theoretical interests iH 0, H, 

aaaaiSEEiiiiHHEiEEa. Despite 

the pairing mechanism being still controversial, a lot of theo- 
retical works indicate that the orbital degeneracy may play a 
key role in these new family of high-Tc superconductors. Dif- 
ferent from the conventional cuperate superconductors where 
the undoped compound is a Mott insulator, the pure LaOFeAs 
compound is a poor metal. This implies that the orbital de- 
generacy may dramatically affect the properties of the normal 
state. It is well known that one of the most basic model to 
understand the strongly correlated systems is the well-known 
single-band Hubbard model. Although thousands of works 
have been focused on such a deceptively simple model, the 
physical properties are not fully understood except the one- 
dimensional (ID) case where the exact solutions are available 
ifTgll . The exact results indicate that the Hubbard model with a 
filling factor one is a Mott insulator at the zero temperature. It 
is quite interesting to ask whether the ID Hubbard model with 
the same filling factor but with an additional orbital degree of 
freedom is still a Mott insulator in the zero temperature? Aim- 
ing to answer this problem, we study the conductivity of an 
extended ID Hubbard model with the orbital degree of free- 
dom in the scheme of the Bethe-ansatz solution. Unlike the 
single-band Hubbard model where the conductivity is found 
to be zero for any nonzero repulsive interactions, the Hubbard 
model with the orbital degree is found to be a conductor when 
the repulsive interaction is smaller than a critical value and a 
phase transition from metal to an insulator occurs when the 
on-site U is larger than the critical value. 

AID electronic system with the orbital degree of freedom 
can be modeled by 



H = -f 2 + h.c.) + uYj ni,ani,a'. 



(1) 



i,a<a' 



where / - 1,2, ...,N identify the lattice site, is the total 
particle number, L is the system-size, and a - 1 , 2, 4 labels 



the four states of single site, i.e. {1 t, 1 i,2 t,2 i) with 1,2 
being the indices of orbital and t, i being the indices of spin 
ll20l I21I1 . The internal degree of freedom in the Hamiltonian 
([TJ is specified to spin and orbital in present model. The C^^ 
creates an electron with spin-orbital component a on site /, 
and «,fl := C^^Cia is the corresponding number operator at site 
/. The system ^ is assumed with periodic boundary condi- 
tion and (f) is the magnetic flux piercing the ring. The system 
([TJ is the Hamiltonian for four-component systems, and there 
are various discussions on multi-component Hubbard model 
in one dimension |21, 22, 23.. 24. 251. 

The conductivity of a many-body system generally takes 
the form 



cr((L)) — 2nDc6{<jS) + cr,.(oj). 



(2) 



where Dc is the charge stiffness and cr^(w) is the regular part 
of the conductivity. If Dc is finite, the system is a perfect con- 
ductor; and if Dc is zero but cTriciJ) is finite, the system is a 
normal conductor; while if both of them are zero, the system 
is an insulator. At zero temperature, the transport properties of 
one-dimensional systems depend usually on the charge stiff- 
ness. Kohn showed that the charge stiffness can be computed 
from the ground-state energy E{(p) as ll26ll 



Dc 



2" ' 



(3) 



where <p is the external magnetic flux. It is well-known that 
the magnetic flux (p piercing the system with the periodic 
boundary condition can be gauged out by imposing the twisted 
boundary conditions on the system [27, 28]. Therefore, solv- 
ing the Schrodinger equation in the presence of the magnetic 
flux with periodic boundary condition is equivalent to that in 
the absence of the magnetic flux but with a twisted boundary 
condition for the wavefunctions 



,x,+L,---) = e'V(xi 



(4) 
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We restrict our studies in the case of L - N. If U - 0, the 
electrons do not interact with each other and the Hamiltonian 
can be transformed as 

H = -2f J] cos(^ + 4,IL)Cl^Ck,a. 
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At the ground state, the electrons are arranged below the 
Fermi surface according to the Pauli exclusive principle. Then 
the density of electrons is 



L~ In 



4dk, 



where kf is the Fermi momentum and 

TlN 

For the case ofN = L, the Fermi momentum is kf - n/A. The 
ground-state energy is 



E{c/>/L) 



E(<p/L) 



At r 

— I cos(A; + (f)/L)dk, 
^ J-kp 



(5) 



AtL 

[sin(^f + - sin(-^f + . (6) 

n 



At the U - Q case, the ground-state charge stiffness can be 
obtained analytically as 

2 V2 

D, = — — = 0.9003163. 
n 

Therefore, the system is a perfect conductor. While if the in- 
teractions between the electrons tends to infinity, f/ — > oo, 
each double occupation will cost an infinite energy thus the 
each site favors the single occupation. The system is an insu- 
lator Therefore, a quantum phase transition from a conduct- 
ing phase to an insulating phase should occur between these 
two limiting cases. 

For finite U, the Hamiltonian is quasi-integrable if site oc- 
cupations of more than two electrons are excluded. Physically, 
this is reasonable for the present studied case with filhng fac- 
tor one due to the state with more than two electrons on a 
site is energy unfavorable. Following the standard procedure 
[122. .23. .241 . the energy of the system ([TJ is 



£■(0) = -2f y cos kj. 



(7) 



where the quasi-momentum kj should satisfy following Bethe- 
ansatz equations 



Jk,L 



sin kj — /!/, + irj 



n 



Aa - sin ki + IT] 



Aa - sin ki - irj 

M , 

PI fig- Ab+ IT] 

\ \ iia- Ah- irj 



i sin kj - Ah - irj ' 

pr Ag- At + i2j] pr flc- ylg + l 

[ \ Ag - Ab - ill] fic - Aa-i 



b=\ 
M 



b=\ 
M' 



Aa + in 



M" 



b=\ 
M' 



pr jla - l^c + i^J] pr Vd - l^g + I T] 

i \ Ug- Ur- ilrj i I Vd- Un- . 



M' 



pr Vg - jUfc + _ pr Vfl - jUc + 
where rj - U/At, A,fi and v are the rapidities. 



(8) 




FIG. 1: (color online). The curve of charge stiffness versus the on- 
site couplings for L = N. The cross point on y-axis is 2 y/l/n. 



For ground state (i.e., at zero temperature), the k, A,fi, v are 
real roots of the Bethe ansatz equations (O. Taking the loga- 
rithm of the Bethe-ansatz equations, we get 



J / sin kj - Ag 



2nlj - kjL + 4) + 2'^ tan 

a 

2nJg - 2 y tan 

/ 



tan 



1 j He - 



tan 



2l 



2* = 2Z--'(^)-2Z--(^) 



-1 / V d - Ha 



2.a. = 22:.an-(r^)-22.an-(^), 



(9) 



where [Ij, Jg, Ka,Qa} are quantum numbers. Ij takes integer or 
half-odd integer depending on whether M - 1 is odd or even. 
Jg, Kg and Qg take integer or half-odd integer depending on 
whether N-M-M',M-M'- M" and M' - M" are inte- 
ger or half-odd integers, respectively. If = 4n for n being 
odd integer, the ground state is non-degenerate, and quantum 
number are centerred symmetrily around the zero point. 

We numerically solve the Bethe ansatz equations ^ with 
the finite system-size L. The charge stiffness versus the inter- 
action U is shown in Fig. [T] We see that the charge stiffness 
shows a sharp peak with the increasing system size. In the 
thermodynamic limit, the charge stiffness is expected as a step 
function, which takes a non-zero value at one side and zero at 
another side. The sudden jump point defines the critical Uc of 
the phase transition. For the present model, if U < Uc, the 
system is a metal while if U > Uc, the system is an insulator. 
While for the single band Hubbard model, Lieb and Wu show 



U/4t 



1/lnN 



FIG. 2: (color online). The derivative of the charge stiffness, 
some certain coupling U„„ the derivative has a minimum. 
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FIG. 3: (color online). The curve of the minimum of derivative of 
the stiffness versus the system size A'. The data can be fitted as 
£),'„ = -0.374 - 0.098^. One see that the minimum is divergence if 
the system size tends to infinity. 



FIG. 4: (color online), (a) The scaling behavior of the U,„. The 
data of the two-band model can be fitted as U„,/4t = 0.5475 + 
1.19993/ In A'. When the lattice number A' tends to infinity, U,„ 
becomes U,., the critical value of UJ4t reads 0.548 ± 0.005. (b) 
The same scaling analysis is performed for the single band Hub- 
bard model. The data of the single band model can be fitted as 
U,„/4t = 0.00483 + 1.79084/ InN. The critical value of f/,/4f reads 
0.005 ± 0.006, which covers the exact critical point Uc = of the 
single band Hubbard model. 
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U/4t 



that the Mott-insulator transition only happens at the Uc - 
case. 

The critical Uc can be determined by the derivative of the 
charge stiffness. The derivative of the charge stiffness versus 
the coupling U is shown in Fig. |2l We see that the derivative 
has a minimum at a certain The minimum is decreas- 
ing with the increasing system-size. When the system-size 
tends to infinity, the minimum is divergent, which can be seen 
clearly in Fig. [3] From Fig. [3] the value of the minimum of 
derivative of the charge stiffness versus the system-size can 
be fitted into a straight line. Thus the charge stiffness becomes 
steeper and steeper as the system size increases. 

The critical coupling U,,, with finite system-size versus the 
system-size is shown in Fig. ID The data of U„, and system- 
size A^(= L) can be fitted as U„,/4t = 0.5475 + 1.19993/ In A^. 
When the system-size tends to infinity, U„, becomes Uc- The 



FIG. 5: (color online). Charge stiffness with holes. Here the system 
size L = 100 and the particle numbers A' = 84, 92, 100. The charge 
stiffness has a sharp peak only at the case of L = A. 



critical Uc/4-t reads 0.548 + 0.005. We also perform the same 
scaling analysis for the single-band Hubbard model [Fig. 
lib)] and find that Uc/4t = 0.005 ± 0.006. The difference 
between the two models is clear For the case of N = L, the 
multiband might enhance the critical value Uc/4-t. 

Now we consider the system with some holes. For this case 
the Bethe ansatz solutions are solved by choosing suitable 
quantum numbers {Ij, Jg, Kg, Q„}. The charge stiffness versus 
the coupling U is shown in Fig. |5j where the system-size is set 
to L = 100 and particle numbers are A^ = 84, 92, 100. We see 
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that only at the case of N - 100, the curve has a sharp tran- 
sition; while for both cases = 84 and - 92, the charge 
stiffness always takes a nonzero value. This observation is 
consistent with the fact that the system would be a metal if we 
add some holes. 

In conclusion, starting from the Bethe ansatz solutions of 
ID two-band Hubbard model, we study the multiband effects 
on the conductivity. We find that multiband would enhance 
the critical value Uc of on-site interactions of the transition 
from a metal to an insulator, while the critical Uc for the sin- 
gle band Hubbard model is zero. The finite system-size would 
have a 1/lnL correction to the actual value. The orbital de- 



gree of freedom might play an essential role in the properties 
of the electronic systems. These results are helpful to under- 
stand why undoped High- T^c superconductors are usually insu- 
lators, while recently discovered iron-based superconductors 
are metal even without doping. 
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